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The Approximate Self-Consistent Field Molecular Orbital Theory 
(ASCFT) with complete neglect of differential overlap (CMDO) has been 
applied to the calculation of electronic structure for several selected 
molecules in this text. The theory leads to calculated equilibrium 
geometry configurations and the barrier to internal rotation of ethane, 
inversion barrier of ammonia and the barrier for cis to trans conversion 
for difluorodiazine. These results are in reasonable agreement with 
experimental values in most cases. The stretching force constants are 
higher than experimental values by a factor of approximately two. The 
orbital energy levels for ammonia and for oxygen difluoride as a function 
of bond angles were in reasonable agreement with the full LCAO-SCF 
Calculations, A variation of parameters to find a new set of parameters- 
was attempted with vinsuccessful results. 

The Extended Hiickel Theory (EHT) can be applied to aliphatic 
hydrocarbons but its application to non-hydrocarbon molecules is quite 



limited 



LIBRARY 

NAVAL POSTGRADUATE SCHOOL 
MONTGF.EY. CALIF. 93940 



TABLE OF CONTENTS 

Section Page 

I Introduction 9 

II Hiickel molecular orbital theory 9 

III Extended HUckel theory 15 

IV Approximate self-consistent field 17 

Molecular orbital theory 

19 

1, SCF equation 

21 

2. Theory with complete neglect of differential 
overlap (CNDO) 

3* Choice of Parameters of CNDOl and CND02. 

28 

V Discussion of Results 

29 

1. Comparison of the equilibrium configurations 
and force constants, 

2. Water molecule (H^O), 

3. Ammonia molecule (NH^), 

U. Ethane molecule (C 2 H^). 35 

5* Methane and Methylfluoride (CH|^ and CH^F). 38 

6 . Difluorodiazine (N^F^). 39 

7. Oxygen difluoride. 

8 . I^drogen Peroxide, 

9 . Homonuclear diatomic molecules U6 

VI Conclusions i;8 

5o 

Bibliography 

Appendix 

53 

I Computer program for ASCFT. 

9I1 

II Computer output of H 2 O (Sample), ^ 



3 



t 



LIST OF TABLES 



1 . 

2 . 

3. 

U. 

5 . 

6 . 

7. 

8 . 
9. 

10 . 

11 . 

12 . 

13. 

111 . 

15 . 

16 . 



Slater exponents for valence Shell S and P atomic orbitals. 

Bonding parameters'/^ (in e.v. ). 

"Average" ionization potentials Jyll (e.v.). 

Difference between CNDOl and CND02. 

Matrix elements from atomic data, LU (CND02). 

Comparison of equilibrium geometry and force constants. 
Calculation of Rom and^HOH for H2O . 

Comparison of the results of calculations for NH^. 

Change in barrier for NH^ w.r.t. minimizing parameters. 

Comparison of Ethane between EHT.and ASCFT calculation. 

Comparison of barrier for Ethane. 

Comparison of r, eAr and "^©between CH|^ and CH^F. 

Results of experiments and calculations for N 2 F 2 . 

Equilibrium configuration between cis and trans for N 2 F 2 molecule. 

Comparison of barrier between C2ir^O-^} and for N 2 F 2 . 

Equilibrium configviratians and some partial derivatives for 
homonuclean diatomic molecules. 



5 



LIST OF ILLUSTRATIONS 



1 . Orbital Energy of NH3 as a function of the p3rramid height (ASCFT 
and full LCAOSCFT). 

2 . Orbital energy of NH^ as a function of the Pyramid height (EHT). 

3. Barrier for N2F2 as a fxmction of 

U, Orbital energy of [>20 as a function of ZFOF* (ASCFT and full 
LCAOSCFT). 

5, structural parameters of H2O2. 



7 



ACKNOTLEGEMENT 



The author wishes to express his appreciation to Dr. William M, 
Tolies for his continuous guidance and assistance in this project and 
also wishes to thank members of the computer facility of the United 
States Naval Postgraduate School for providing sufficient computation 
time. 



8 



I. INTRODUCTION 



There are many levels of sophistication for the calculation of 
electronic structure using molecular orbital methods. The simplest 
approach is the Hiickel Molecular Orbital Theory (MO). This is one of 
the semiempirical methods which is used extensively in the quantum 
mechanical study of the mobile electrons of conjugated molecules. The 
HMO was extended by many people. One of the theories examined in this 
paper is the '‘Extended Hiickel Theory (EHT)" which was developed by 
Roald Hoffmann (1). A higher level of sophistication is the self- 
consistent field (SCF) molecular orbital calc^^lation, One of those 
methods which has recently been developed is based on a theory intro- 
duced by J, A. Pople: "Approximate Self-Consistent Field Molecular 

Orbital Theory (ASCFT)" (2), (3), ih). This ASCFT Calculation is 
predominantly used throughout this text. 

We can say molecular orbital methods are "curve fitting". They 
are really based on semiempirical methods. The purposes of the calculation 
described herein are the prediction of potential surface by Pople 's 
theory to examine and to try to improve the agreement of a calculated 
potential surface with an experimental potential surface. 

The background describing each of these theories will be discussed 
in the next section. 

II, Httckel Molecular Orbital Theory (Ht-IO) 

In this theory, the mobile electrons can be treated as occupying 
a set of delocalized molecular orbitals (by the Pauli Exclusion Principle, 
not more than two electrons in each molecular orbital). These molecular 



orbitals are eigenfunctions of a one electron Hamiltonian, The 
orbitals are approximated by a linear combination of atomic orbitals 
(LGAO) centered on various atoms. The Slater type atomic orbitals 
(^) are generally used as approximations to the actual atomic orbitals 
for each atom. When this theory is applied to a wide range of molecules, 
it has the advantage of great simplicity but serious discrepancies result 
because of the neglect of interelectronic repulsions for the case of 
many electrons. 

For the unsaturated molecules : The JC - molecular orbitals are 

formed by a LGAO of the 2ft atomic orbitals which, it is assumed, share 
the same nodal plane. In hl“X electron systems, the molecular orbitals 
are given by 



( 1 . 1 ) 



uf, = iCf- 4 

where represents the molecular orbital, ^^^represents the atomic 

orbital in the 2^ atom and represents the coefficients of the 

fh 

molecular orbital and the atom (atomic orbital). It is noted that 

the subscript i, j, k, 1 etc, stand for the molecular orbitals and 
^ A etc, stand for the atoms or atomic orbitals. 

From a given electron Hamiltonian, Equation (l,l) must satisfy 
the equation 



m-, -- 



( 1 . 2 ) 



Where Ej is the eigenvalue of the j 
and'^ is the Hamiltonian operator. 



th 



state, is a column matrix. 
This can be put in matrix form; 
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( 1 . 2 )' 



We multiply equation (1.2) by \Jj^ and integrate over all space, dr , 
then we obtain 






(1.3) 



Substituting equation (l.l) into (1.3) 



yua/ 

Defining 

and 



3/ ' /q O' 



Hu 2. = f <k, dr 



S/j» dr 

¥e may write equation (l.li) 

^ Ca/i (^p - E Cp CpS/ii/ 

Taking the derivative of equation (1.7) with respect to Gfc\ to 
the energy(^ f obtain 

It may be shown that 

~ Si;yU 

H/43; - 

Therefore we obtain the equation 

51 Ci^ IjqA = 5^^ 



(1.5) 

( 1 . 6 ) 



(1.7) 
miniinize 

( 1 . 8 ) 

(1.9) 

(1.10) 



( 1 . 11 ) 



or 






( 1 . 12 ) 
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This equation (1.12) is the well-known secular equation and gives rise 
to the secular determinants 




(1.13) 



Equation (1.13) is the non-trivial solution of (1.12). 

Now we discuss the Hermitian Property of the Matrices ( 
and A matrix H is said to be Hermitian if it is equal to the 

conjugate of its transpose, i.e., 

H' H 

Whens H = ( 

( I; 3 0 

By definition, a real Hermitian matrix is symmetric, but an imaginary 



(I.IU) 



symmetric matrix is not necessarily an Hermitian matrix. It is noted 
that the eigenvalues of an Hermitian matrix are all real. 

Next we discuss some transformation properties. In general, the 
atomic orbital does not form a complete orthonormal set. This 

basis, , can be transfonned into an orthonormal basis by using a 
unitary transformation matrix such that 

■ new 

^ (1.15) 

where ^ 2 j the new orthonormal basis formed by transforming the old 
basis 4)pn in N dimensional space and is the transformation matrix 

element. The old H can also be transformed into a new basis such that 

hW - Z j (1.16) 

It may be recalled that the old basis is 



Therefore equation (l.l6) can be rewritten 

,Tie\0 t — 7 ■* . , old- 



^ ^ . I oia-p 

•>323 " Zj \a<7 Ho-A 



(1.16)' 
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or 



new ^ ★ old / 

r1ju2/ - ^ •Aio‘ \>z/ 



It is noted that if the transformations were made: 



r = ? 



(1.17) 



(1.18) 



then the equation (1,17) becomes 

h;: = z, 

^ (Ti?\ ^ 

Equations (1,17) and (1,19) can be written by matrix notation such that 



(1.19) 



and 






(1.17)' 

(1.19)' 
then we have 



From equation (1,19)', if T is a unitary matrix^”}” — 

|^new_ y-i 

This equation (1,20) is frequently used and usually T is real and 



( 1 . 20 ) 



is a symmetric matrix. 

Referring to the secular equation (1,12), we define^-j^l^s^, then 

I ( = 0 

This can be expressed in matrix notation in the following manner 



HC-SC/\-0 

;^=art) 

= (Cm) 

s-.(s>/i (;>>,/( ,i- 1 -—m) 

Therefore, we obtain 
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It is generally true that the matrix (S”^H) may not be symmetric. In 
equation (1.22), the eigenvalue cannot be obtained by methods usually 
employed to diagonalize Hermitian matrices. One way to solve this 
problem is to transform the matrix S to an orthonormal basis, then 
diagonalize the H matrix. The procedures are described as follows: 
Using the transformation matrix Q such that we transform the 

matrix S to the identity matrix I 



5"“- (q'qV 



Also we diagonalize the matrix S by a unitary transformation matrix T 
such that 

> = TS T' (1.25) 

or - (T'X'TT' 



-I 



( 1 . 26 ) 



= [(T^'^Ty(T'>i''^T)] 

Equating equation (l.2ij.) and (1,26), we obtain 

0=T^/)’^T 

Then we transform the old basis H to the new basis by the transformation 
matrix Q such that 






(1.28) 



Finally we can diagonalize by the matrix U, Then we obtain the 

eigenvalues 









(1.29) 



lU 



where U is the matrix which diagonalizes the and T is the unitary 

matrix which diagonalizes S, Also the correct eigenfunctions are 



The general operations for the application of HMO theory to 
specific problems are 



III. Extended Hiickel Theory (EHT) 

Hiickel calculations originally were applied exclusively for the 
Planar conjugated and aromatic systems of JC electrons with the neglect 
of overlap for non-ad jacent atoms. In the extended Hiickel theory (EHT), 
Roald Hoffmann (1) expanded the basis set to consist of 2s and 2P 
orbitals for Ceirbon atoms and Is orbitals for Hydrogen, with the 
inclusion of all overlap integrals. The EPIT can calculate molecular 
orbitals for the aliphatic and aromatic compoiinds which gives good 
results for most hydrocarbons because of the res enable set of parameters 
chosen for the calcvilations , The parameters are nicely adjusted to the 
self-consistent value for the hydrocarbons but when the EHT is applied 
to non-hydrocarbons, it is not as successful, and even predicts a lack 
of minima in many bond stretching motions. 

III-A. Method of Calculation 

For the calculation of the molecule CnHm, we use a basis set 




(1.30) 



(a) Setting up the Secular determinant 



(b) Expanding the Secular determinant 



(c) Finding the eigenvalues 

(d) Determining the M.O, *s by finding the eigenvectors. 
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consisting of m-Hydrogen Slater orbitals, n-2s and 3n-2P orbitals. 

The Slater exponents v;ere selected to be 1,0 for Hydrogen and 1,625 
for Carbon throughout the calculations. It is noted that 1,20 was 
selected for Hydrogen in the calculation of the ethane raolecule. Those 
values are listed in Table 1, The order of the matrix is Un + m which 
gives the number of energy levels of the system and the eigenvalue E's 
for off diagonal matrix elements are retained in the Secular Equation, 







TABLE 1 






Slater exponents for 


valence Shell S and P 


atomic orbitals. 


H 


1.0* 


C" 


1.U5 


0" 


2.10 


Li 


0.65 


C 


1.625 


0 


2.275 


Be 


0.975 


C+ 


1.80 


0“^ 


2.1;5 


B" 


1.125 


N 


1.95 


F 


2.60 


B 


1.30 




2.125 






* Present 


calculations 


we carried out using 


1.2 


for Hydrogen. 



The diagonal elements are chosen to be valence state ionization 

potentials 3 the particular values used are essentially those of Skinner 
and Pritchard (6), For example, the SP^ valence state ionization 
potentials for Carbon and Is Hydrogen atoms have the following values 

H/j/u (C2p) c — 1 1 . 4 

(Gs) z -2\< A ( 2 , 1 ) 




The choice of the off-diagonal term 




e.v. 

is approximated 



as 
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T"]aI2/ - K ( ^/yA + H 3 / 2 ) (,SyA3^ 

The approximation (2,2) was first discussed b 7 Miilliken (7) and used 
in the M.O, Calculations by Xfolfsberg and Helmholtz ( 8 ). 

The value of K in equation (2,2) was chosen to be 1,75. These 
calculations give gairly good results for hydrocarbons. 

IV, Approximate Self-Consistent Field Molecular Orbital Theory (ASCFT) 

Approximate methods for obtaining self-consistent field moleciilar 
orbitals have recently been developed by J, A. Pople (2), ( 3 ), (U), 

Slater type orbitals (5) are used as basis functions. Overlap integrals 
between two different atoms were calculated using the formula given by 
I'fiilliken, Rieke, Orloff and Orloff (9), It is assumed that the Slater 
type atomic orbitals are real atomic orbitals, and that the M.O. 's for 
the groimd state are obtained by Linear Combinations of Slater type 
atomic orbitals. If the Coefficients in the LCAO Molecular orbitals 
are chosen to minimize the total energy, one obtains LCAO - SCF orbitals. 
The best LCAO approximation to SCF functions were first given by Roothaan 
(10), It is noted that the EHT treated all the electrons independently, 
EHT wave functions are used as a starting point for the self-consistent 
field interations in this ASCFT, 

The SCF methods were developed using two approximations based on 
the neglect of the differential overlap between all valence atomic 
orbitals. One of these involves complete neglect of differential overlap 
(CNDO) in all basis sets connected by the transformations as mentioned 
in Section I, that is: the off-diagonal elements of overlap matrix 

in the Hamiltonicin matrix are to be neglected. The other involves the 
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neglect of the products of of any two atomic orbitals 

and in all electron-repulsion integrals. This complete 

neglect of differential overlap of all basis sets will result in non- 
orthonormal basis sets. 

The CNDO method in ASCFT should be regai?ded as an approximation 
to a full LCAO - SCF calculation which uses a minimal basis containing 
only inner-shell and valence-shell atomic orbitals. The Slater 
exponents used in these calculations are listed in Table 1 (for Ifydrogen, 
1,2 is used). The values for ionization potentials in 2s and 
2p atomic orbitals from Lithium to Fluorine were set semiempirically. 
Those values used in these calculations are tabulated in Table 2 and 
Table 3 respectively. 



TABLE 2 



Bonding Parameters /u a (in e.v. ) 



Atom 



H 

9 



Li 

9 



Be 

13 



B 

17 



C 

21 



N 

25 



0 

31 



F 

39 



TABLE 3 

"Average” Ionization Potentials (in. e.v.) 



Atom 


H 


Li 


Be 


B 


C 


N 


0 


F 


Is 


13.06 




— 


— 


— 


— 


— 


— 


2s 


— 


5.39 


9.32 


lU.05 


19.UU 


25.58 


32.38 


U0.20 


2P 


— 


3.5U 


5.96 


8.30 


10.67 


13.19 


15.85 


18.66 
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The CNDO method of ASCFT can be used for closed ground shell 
structures and doublets, but a triplet structure is not available with 
the computer program used (for example P 2 , ^ and C 2 ). The computer 
program is presented in Appendix 1. 

IV-1. Self-Consistent Field (SCF) Equation 

The initial discussion covers Roothann's LCAO - SCF method (10) 
as a basis for calculating the M0‘s of a molecule, paying particular 
attention to transformation properties. The Is Hydrogen and valence 
shell electrons for Li to F are assigned to LCAO molecular orbitals. 
Minimization of total energy by a variation principle with SCF orbitals 
leads to a series of simulteineous non-Linear equations, the so called 
Hartree Fock equations (11), which are applied to LCAO molecular 
orbitals with a closed shell configuration. These equations were shown 
by Roothaan to reduce to the form (10) 

where (using atomic units) 



^2^ = + (j/iZf 


(h.2) 


II 


(U.3) 




(h.h) 


r| [^(0 -pL 


ih.$) 
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occ . 



(U.6) 

(U.7) 



Where Ei is the orbital energy for the i^^ M.O. H/JU^ is the matrix 
element of the one electron Hamiltonian including the kinetic energy 
and the Potential energy in the electrostatic field of the core, this 
being written as a sum of potentials V^(r) for various atoms A in the 
molecule, is the matrix element of the potential due to valence 

electrons and depends on the M.O. ’s via the population matrix R\cr . 

The total energy of the molec;ile is obtained by adding the repulsion 
energy between isolated cores and the separated valence electrons. The 
repulsion term can be approximated by a point charge model 



Etotal = 

A<e> ^ 



(U.8) 



where is the core charge of atom A. The electronic energy is 

given by 



-Ka|) 

-- L H/.Z, + -^2 



(U.9) 
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(li.io) 



Select = 2 H/«2< + 

A2^ 

+^ 2 >] 






17-2* Theory with Gcanplete Neglect of Differential Overlap (CNDO) 

First consider the simplest version of the theory which involves 
the complete neglect of differential overlap between atomic orbitals 
on the same atom. The approximations are given below. 



The atomic orbitals are treated as if they form an ortho- 
normal set such that 3||2; “ then the 0^ forms an 
orthonormal matrix. The condition for orthonormality becomes ; 



Approximation A 



A CviQ-m = 



i 



^ ^ Q Ol.ll) 

With this approximation, the diagonal matrix elements (electron 

density) correspond to the electron populations of the atomic orbitals 

4Uj 

l^yU = (U.12) 



where 2N is the number of valence electrons. 

Approximation B . 

Equations (li.it) are neglected because of the orthonormality of 
atomic orbitals. This means ( ) is zero unless and 

. The non-zero values will be written such that 
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(U.13) 






Appraximation C 

The electron- interaction integrals are assumed to depend 

only on the atoms to which the orbitals and belong and not 

on the actual type of the orbitals. This means that there remains only 
a set of atomic electron-interaction integrals measuring an 

average repulsion between an electron in valence atomic orbital on atom 
A and another in a valence orbital on atom B, The set of atomic orbitals 
on one particular atom A is replaced by an alternative set tm. 






~ ^ Owyu 

y (li.lU) 

A general electron-re plus ion integral in this new basis set is given by 



Now we apply approximation C, then equation (U.l^) becomes 



z (h.l6) 

It is noted that equation (l|.l6) is invariant under the orthogonal 
transformation matrix 0 which gives the new basis "t. 

Using approximation A, and B, the matrix elements become 




H"/* ^ 

Ham 



(li.l7) 
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Zf 



(U.18) 



Equation (I|.»17) and (U»l 8 ) can be rewritten by approximation C 

(m) 



\^-2f ~~ ““ O 



0i.l9) 



(ii. 20 ) 



where yM belongs to atom A and P 33 is the total valence electron 
density on atom B, 



Pbb = E 

Also the core matrix ^ separated as follows: 

H*>u - 



(U.21) 






BCfA) 



(U.22) 



B(#4) 

there \J/iyi is the diagonal element of 4v with respect to the 
one-electron Hamiltonian containing only 

the core of its own atom ( U/j/j is an atomic quantity measuring the 
energy of the atomic orbital) and is the one-electron matrix 

element using the local core Hamiltonian ( JU and 29 are on the 
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same atoms). It is noted that the cores consist of the nucleus and 
inner-she]J. electron. 



Approximation D 

The interaction terms of the distribution on atom A 

and the nucleus of atom B, ( ) are equal to zero if /Jl 

The integrals (J>\Va\u ) are taken to be the same for all valence 
atomic orbitals on Atom A and there results 






(U.2U) 

It should be noted that the attraction term matrix element V^b 
between nucleus B and the adjacent valence electron ytj in the field 
of the nucleus A need not be symmetric. Therefore the equations (k*22) 
and (ii.,23) can be rewritten 



- \J//^ ( ja 



(U.25) 



~ 0 



Jflhic. j 



(ii.26) 



It is noted that if and 4^^ are on different 

atoms which will be discussed in the next approximation. 



Approximation E 

If and are on different atoms, it is only necessary to 

consider the possible lowering of energy levels by being in the electro- 
static field of two atoms simultaneously which is referred to as a 
'•Resonance Integral" and is denoted by can be 



2h 



approximated as 



~ 



(ii.27) 

where ySfiZ is a parameter depending only on the nature of the atoms 
A and B. Empirical values will be used for these parameters. 

The f^: 2 / matrix elements now reduce to the form from equation 

(h.l 9 ) to (I1.27) 

^ B(^A) ^ (U.28) 



(I1.29) 



Ku2) ^ 2 ; — l^2;^3 

I 1 ( 

where (|^ belongs to atom A and belongs to atom B. 

Using the same approximations, the total energy is given by one 
or two atom types from the equation (U»8) 



Cw.l -LEa ZZ Eab 

^ A<e> 



where 



AA 



(U.30) 



Ea 6 ?/W^aZ) 

BiaPbb'^^B ""RaI^B“Pb6\4>1 



0i.3i) 
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rV-3. Choice of Parameters of CNDOl and Cl® 02 

In the complete neglect of differential overlap approximation, 
the CMDOl and CND02 methods are different levels of approximation 
of the equations. 

The LCAOSCF Hatree-Fock matrix elements are 




(U.28) 



(h,29) 



where the bonding parameter is given by 

In the equation (U*32), i® chosen empirically but depends only 

on the nature of atcM A, These values are listed in Table 2, The 
parameters Li?5j25 and U^pj2p can be estimated from ionization 
potentials, using the theoretical value of IfxX (12). 

Is(X; . £{X]2s'if) -£(X, 2s''2f) 

= -U2S,2S - 

Ip(X,2s"'2p'’) = £(Xt2s>”">E(X,2S"2p") 

= — U2p^2P — ("Wfh-0 

or 

lyM = (U.35) 

and the atomic electron affinities AyU S are used by the formula 



(U.3li) 



■AyU = U/yA -^'Za'^AA 



(lt.36) 



26 



W|H W|H 



Also the equation (U.28) can be rewritten 



B(M) 



(U.28)' 



where the last terra of equation (U.28)* is described as penetration 
integral contributions to 



The following table gives the approximate form between CfJBOl and 
CMD02. 

TABI£ k 

Difference between CMDOl and CND02 



reference 


CNDOl 


CND02 


XJ/AJLA 


lyM +(Za-|)^A 


“g" 4 ^ 






(U..37) 






U/J/i 


¥here 1 


+ Z(^e'!!4B“ ^b) 

im ^ ^ '<=*''1 (U.36) 

the average ionization potentials (e.v.) used to fix 

\Jujji are listed in Table 3* 


2 


the matrix elements §(ls + As) 
in Table 5. 


Sind J(lp + Ap) are listed 



TABLE 5 

Matrix Elements From Atomic Data (e.v, ) 





H Li 


Be 


B 


C 


N 


0 


F 


^(la+Aa) 


7.176 3.106 


5.9U6 


9.59U 


111. 051 


19.316 


25.390 


32.272 1 


Klp+Ap) 


1.258 


2.563 


U.OOl 


5.572 


7.275 


9.111 


11.080 
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If the penetration integral term in equation (U»28)' was neglected 
to avoid the excess bonding and high value of force constants, then 
equation (1|,28)' can be reduced to the form of (U.38). This corresponds 
to substituting the following term for CND02 from CNDOlt 

Using parameters specified in the above way, the LCAOSCF equations are 
solved in the following steps j 



(a) An initial set of M.O, coefficients CjCyti is obtained from 
a'Huckel type” calculation such that is replaced by 

average ionization potential and • 



(b) 


Electrons are assigned 


in pairs to the M,0, 's with 


lowest energies 


(c) 


Calculated and ] 


are used to form a new Hatree Fock 




matrix element 


• 




(d) 


New coefficients Q^is 


are derived from this 


and the 



process is repeated from step (b) until self-consistency is 
achieved on all coefficients with a tolerance of 0,0001, 

V, Discussion of Results 

The molecviles selected for these calculations are as follows: 
Hydrogen molecule (H 2 ), Lithium molecule (Liz), Beryllivtm (Bez), 
Nitrogen molecule (Nz), Fluorine molecule (Fz), water molecule (H 2 O), 
Ammonia molecule (NH^), Ethane molecule (C 2 Hg), oxygen difluoride 
(F 2 O), Difluorodiazine (N 2 F 2 ), Methane (CH|^), Methylfluoride (CH^F) 
and Ifydrogen Peroxide (H20^), These calculations will show the minimum 
configurations for bond lengths, bond angles, force constants, the 
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energy level ccaisiderations for NH^ and F2O, and the Barrier for 
C2^^f N2F2 and H2O2. Finally the calculations will extend the variation 
of parameters to fit the experimental potential surfaces. 

V-1, Comparison of the equilibrium configurations and force constants 
The ASCFT calculations with various inter-nuclear distances and 
bond angles are useful in predicting equilibrium bond lengths and bond 
angles. These results are listed in Table 6. The equilibrium configur- 
ations (bond lengths and bond angles) are in fairly good agreement with 
experimental values but the force constants are too high compared with 
the experimental values by a factor of about 2.0. This can be improved 
by adjusting parameters, such as the Slater exponent and yO . 

Sample results are discussed in Section V-9. 



TABLE 6 

Comparison of Bond Lengths, Bond Angles & Force Constants ("^r ^ 




Bond Length 


('l^'lForce \0^ N 

Bond Angle (deg) '' -^'Consts V-k©: 1(5" erg/md? ) 




Cal. Obs . 


Cai. Ubs. 


Cal. 


Obs. 


H2 

N2 


0M9 0.7U15 




10.11i3 


5.6 

/’n), \ 


1.1395 1.0976 




52.160 


(.-in/ 

23.10 

- - ,, - . - ... - .» 


H2O 


1.0289 0.958 


10U.22 10U.U5 


17.260 




0.8U5 


0.69li I 


NH3 


1.0657 1.015 




lh.68 






106.825 106.6 


0.95 




0% 


1.1135 1.1135 




12.65 


5.ii5 




109.U7 109.U7 






Rch 


1.1195 1.102 




12.12 


5.35 


I.U633 1.5U3 




l6.3li 


U.57 . - 




157775 io9?37’ 




0.51 


j 


110.86 109P37’ 




0.63 
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V-2. Water molecules (H2O) 

For H2O, C21T symmetry was assumed with an 0-H bond length of 
0.96A°. The first calculation shows an energy minimum at a bond angle 
of 107»28°, At this bond angle, the equilibrium 0-H bond length of 
1,0289A° is obtained. With this bond length, the second calculation 
shows the bond angle of 10ii,22*^ at VThich the 0-H bond length is 
calculated to be 1,0289A°. This is in satisfactory agreement with the 
experimental bond angle of 10U,U5° (l3)» These calculations demonstrate 
the 0-H bond stretching and angle HOH bending motion such that the 0-H 
bond lengths tend to have symmetrical stretching motion and the bond 
angle tends to have scissoring motion (symmetrical in-plane deformation). 
A series of these calculations is listed in Table 7. The stretching 
force constant is calculated to be 17.26 x 10^ (*^^®/cm) which is 
higher by a factor of 2,1 with respect to the experimental value of 
8,23 X 10^ C^y^®/cm), The bending force constant is calculated to be 
0,8U5 X 10“^^ (®^S/rad2) which is higher by a factor of 1,22 than the 
experimental value of 0,69U x 10“^^ (®^S/rad2) ( il| ), 

TABLE 7 

Calculation of Rqh and HOH for H2O 

Calculation I Calculation II Experiment 

Roh(A°) 0.96 (Exp) 1.0280 0.9?8 

zHOH 107.28 10U.22 10 U.U 5 

* The computer outputs (sample) are listed in Appendix 2, 



V-3. Ammonia Molecule (NH3) 

The ammonia molecule with an NH bond length of l.O^A® (13) is 
calculated to be non-planar with an Z. HNH bond angle of 106,82^° 
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which is in close agreement with the experimental value of 106.6° (13)« 

At this value of the bond angle, the NH bond length was fo\md to be 
1,0657 A°. The subsequent calculation shows a bond angle of 10^,86° 
and again the bond length of 1,0673 A°, On the other hand the equili- 
brium bond length for the planar configviration gives 1.05U3 A°, It is 
interesting when calculations presented here are compared with results 
recently calculated by R. Moccia and L, Randaccio (17). They reported 
the minimum bond length of 1.878 a.u. (0,9938A°) at a bond angle of 
105°36' and a bond length of 1,837 a.u, (0.9721 A°) at a planar 
configuration. The change in minimum bond length between non-planar 
and planar configuration is 0,0217 A° by Moccia and O.Ollli A° by the 
ASCFT, The smaller value of this bond length change in the present 
calculation is due to the higher stretching force constants for ASCFT, 
Also, subsequent calculations can predict the symmetrical bond stretching 
motion and symmetrical in-plane deformation of bond angle for the ground 
state of this molecule. The stretching force constants are calculated 
to be lii«6li X 10^ (^^°/cm). This differs by a factor of 2,lh with 
respect to the experimental value of 6,86 x 10^ (lli) and the 

bending force constants are calculated to be 0,95 x 10"^^ erg/rad. ^ 

(factor of 1,58 with respect to experimental values of 0,6 x 10~^ 
erg/rad^ (llj.)). 

It is also interesting to compare the barriers calculated by the 
several methods and the difference in barriers with respect to adjusting 
the various parameters. These are listed in Tables 8 and 9, respectively. 
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TABLE 8 

Comparison of the results of Calculations for NH3 



Reference 


Method of calculations 


Rnh fdU) 


^ HHH (J^ 


Pyramid Barrier 

Height(a,u. )Height(a,u,) 


Present Calc. 


Approximate (CNDO) 
SCF MO 


1.9923 


106 °U9’ 


0.7U9 


0.0251 


Higuchi (18) 


Approximate polycenter 
SCF MO 


1.916 


108.0 


0.592 


0.01U7 


Kaplan (19) 


Single Determinant 
polycenter SCF MO 


1.916 


106°1;7, 


0.719 


0.1ii6 


Moccia (20) 


Single Determinant 
one center SCF MO 


1.928 


110°2t 


0.62U5 


0.0027 


Joshi (21) 


Single determinant 
one center SCF MO 


1.867 


109°3U> 


0.626 


O.Oiill 


Moccia (17) 




1.878 


105°36» 


0.739 


0,00905 


Experiment 




1.9162 


106 °U7 ’ 


0.719 


0.0092(^2) 



From Table 9 , it is interesting to compare the change in barrier 



with respect to a change in bond angles at the fixed bond length of 
1.06$7A° (Calculation II and IV) such that the change in barrier is 
0.000030279 a.u. (O.OO25U083O - 0.002510551) and the change in bond 
angle is 1.965° (106.825° - 10U.86°), From this calculation, it is 
found that the calculated barrier is O.OO25U083O a.u. which is different 
by only 1.2^ from the first calculated barrier. 

The electronic configuration is found to be 

(laj^)2 (le)^ (2a^)2 ^ ^ 

Configurations calculated by ASCFT, EHT and full LCAOSCFT by H. D. 

Joshi (21) are compared in Fig, 1 and Fig 2, There is complete agreement 
between ASCFT and full LCAOSCFT calculations for the orbital energies 
as a function of height of fyramid. It is noted that the lowest orbital 
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Torn ENEI?&Y((li a) 



T^. 2 . Etol.&£ow..oF NH3 vs. PtWlID HU&HT 
(BY EHT) 




3h 



TOT/»L ENERa-Yfe.v.) 



energy in LCAOSCFT calculation is due to the (is) electron (inner- 
shell) of Nitrogen but ASCFT neglects that part. On the other hand, 
the EHT calculation (Fig. 2) does not even demonstrate a double minimum 
potential curve. Therefore the EHT is quite inferior for this kind of 
molecule. 



TABLE 9 

Change in Barrier for NH3 w.r.t. minimizing the parameter 





RSh (.A°) 




2. HNH 




Bcirrier (a.u^ 


Calculation I 
(General) 


1.02 (fixed) 
1.02 (fixed) 


0 

112.0 


120 

106.825 


-13.85761125 

-13.87775028 


0.02013903 


Calculation II 
(Equilibrium) 


1.05U3 

1.0657 


0 

112.0 


120 

106.825 


-13.86Ii30636 

-13.889Ul2ii5 


0.02510551 


III 


1.05143 

1.05143 


0 

113.18 


120 

105.10 


-13.86I430636 

-13.8889015U 


0.02li535l8 


17 


1.06^7 

1.0657 


ii3!55 


120 

10U.86 


-13.86362706 

-I3.8897IU66 


0.02608760 


Exrorimental 




112.206 


106.6 




0.0092 



« The , true„Barrier for this calc 
is (13. 88971^66 - I3.86C3O636) 



alculation 
(a.u. ) 



0.025ii0830 



V-U, Ethane molecule (C2H^) 

Calculations based on ASCFT and EHT have been carried out for 
Ethane in both eclipsed and staggered configurations with experimental 
values such that = 1.5U3A°, = 1.102A° axi^J-CCW = 109°37« 

(tetrahedral angle). The first calculation shows the equilibrium angle 
for both eclipsed and staggered forms using EHT and ASCFT, The energy 
then was minimized by varying the bond length ^CC^ and finally by varying 
bond length The stretching force constants ( and ) were 

found and those values are listed in Table 10. 
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It is of interest to look at the barrier resulting in the change 
in geometry. These changes in parameters and the corresponding barriers 
are compared in Table 11. These comparisons show that the eclipsed form 
is more open than the staggered form by a few tenths of a degree; also 
the bond length is elongated a few thousandths of an Angstrom and 
the bond length Rgg is shortened a few ten thousandths of an Angstrom, 

Some rewearchers have shown an elongation for R^. Although these changes 

in geometry are small, their effect on the internal-rotation 
barr ier is very significant. 



TABLE 10 

Comparison of Ethane Betn. EHT & ASCFT Calculation 







Calculated 








Reference 


ASCFT 


EHT 


Expt. (16) 




Eclipsed 


Stagger 


Eclipsed 


Stagger 


D3d 


^CCH 


111,28 


110,86 


110.99 


llO.Uli 


109°37 » 




107,60 


108, OU 


107.90 


108. i|8 


109°37' 




l,ii660 


l.ii633 


1.8070 


1.8000 


1.3U3 


Rch 


1.1193 


1.1193 


0.9086 


0,9088 


1.102 


-^zccHrio" ■ 


(1.7U) 


(1.67) 


(8.83) 


(8.8) 


0.63 


•^HCH [ 


( ) 


( ) 


( ) 


( ) 


0.31 


-krc r,o^ ■ 


13.887 


16.3UO 


197.0 


I8ii.0 


U.37 


L . 


12.123 


12.12U 


16U.03 


186.36 


3.33 


3arrr<sr . 


2.3288 


2.6386 


2. 873!. 123 
(23) 



There is some uncertainty in the Table 11 for EHT calculation. 
Every reported calculation scheme has the same tendency to increase 
the barrier when the parameters are changed in sequence, but EHT 
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calculation shows the barrier is decreased. The resultant barrier by 
EHT is 1.6358 ^^^/mole. This is not in agreement with the value 
2.6586 ^®^/mole found by varying only one parameter, ^CCH , holding 
the other two parameters R” and R~ constant. It is noted that 
R, Hoffmann (1) reported U.O ^^®^/mole with a Slater exponent (H) of 
1*0, Present calculations use the value of 1,2. However the results 
are in fairly good agreement with experimental values using both ASCFT 
and EHT (with scane uncertain tj^ 







TABLE 11 










Comparison of Barrier for Ethane 






Reference 


Parameter (Eel . — Stag ) 


Barrier (E c lip, 


Stag) 




^CCH ((Aeq'i 


Kec CA*) 


KCH (A^) 


ZCCH 




Rch 


T.L.Hill^2^5 

(I9li8) 


0,3h 


0.005I1 


0.0002 


0.613 


0.665 


0.700 


L.S. Bartell'^^^^ 
(I960) 


0.20 


0.0029 


0.0001 


0.367 


0.383 


0.39U 


J.T.Vanderslice''^^^ 
& E.A.Mason(l960) 


0.17 


0.002U 


-0.0001 


0.387 


O.U16 


0.U3U 


L.Salem(27 ) 

(1962) 


0.05 


0.0007 


0.0000 


o.ioU 


0.106 


0.107 


K.E.Hovilett^^®^ 

(i960) 


o.Uo 


0.0059 


-0.0002 


0.779 


1.108 


1.298 


(29) 

V.Magnasco^ ' 
(1962) 


0.55 


0.0080 


-0.0011 


-1.532 


1.675 


2.9Ui 


0,J, Severs & 
M.Karplus (1966 ) 


0.20 


0.0029 




2.311 


3.611 




Present 

(ASCFT) 


0.ii2 


0.0027 


-0.0002 


1.8637 


2.28U8 


2.3288 


(EHT) 


0.53 


0.0070 


-0.0002 


2.6586 


1.57ii6 


1.6358 


Expt. 3 


.03±0.3(^°^V"iole) (31) 


2.875!0.125 (^°^Vmole) 














(23) 
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Methane (CH^^) and Methylfluoride (CH^F) 

Calculations between methane and methylfluoride were carried out 
using tetrahedral symmetry for methane and C3v symmetry for methyl- 
fluoride with R^h = 1.1135a® for C% and R^ = 1.109 and R~ = 1 , 38 ^ 
for CH^F. The calculated values for the equilibrium configurations are 
listed in Table 12. The force constants are calculated for each molecule 
and these are listed in Table 12. The stretching force constants are 
different by a factor of 2.3 for CH|^ and 2.5 for CH^F with respect to 
their observed values (l5). 

For a methane molecule, the calculated bond lengths and bond angles 
are the same as the observed values. 

For methylfluoride molecules, the bond angle Z HCH is more open 
than the angle of ^ HCF which can be explained as follows s there exists 
rehybridization for the CH^F molecule because there miist be more 
P-character in C-F bond due to the electron withdrawing fluorine; on 
the other hand, there must be more S-character in C-H bond which gives 
more SP^ hybrid character for the C-H sigma ( CT) bonds. It is obvious 
that the SP^ hybrid orbital is more planar than SP3, Therefore <^HCH 
is more open. It is here that the best agreement between theoretical 
and calcxilated results was found. 
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TABLE 12 



Comparison of T > 0 and kr ^ k© between CHj^ & CH^F 





Bond Length (A°) 


Bond Angle (deg) 


Force Consts, 




Cal. 


Obs. 


Cal, 


Obs. 


Cal. 


Obs, 


CH^ 


1.113^ 


1.1135 






12.65 


5.U5 






109. U7 


109. U7 


(.7U5) 


0.1;98 




1,U80 


1.109 






12.09 


U. 71(15) 




1.3U58 


1.385 






23.57 


5.6Ii(l5) 


<'^3':hch 






110.382 


109.997 




0.UU3(15) 


^cri 






108, 5U 


108, 9U 


(1.77) 


0,792(15) 



V-6. Difluorodiazine (N2F2) 

Calculations have been carried out for difluorodiazine in both 
trans and Cis configurations with Rnn = 1,2^A°, = l,l;i^A° and 

= 11^® (13) • In both configurations, the new equilibrium 
distances and bond angles are determined and are listed in Table II4.. 
When canparisons are made of the equilibrium geometry in Table lU, the 
bond angle is more open for Cis than trans configurations. 

Significantly, the bond length Rjjjj is elongated a few thousanths more 
for trans than Cis (same for RmF ). 

It is of interest to compare the energy level difference between 
Cis and trans configurations. The experimental results giving the 
stability between both configurations are listed in Table 13, In 
1963 > Armstrong and Marantz published that the Cis form is more stable 
by 3,0 i 0,3 (36) (refer also (37)), The present calculation 

predicts the Cis form to be more stable by 3»2U9 *^^^Vniole than the 
trans (C2h) form. Fortunately, these results are in most satisfactory 
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agreement with experimental value. In Table lii are compared the diff- 
erences in energies between the Cis and trans configurations for a 
variation in molecular parameters. First the angle .ZNNp was varied 
to minimize the energy. With this bond angle, the minimum bond length 
was calculated and finally the minimum bond length Rup was found. In 
each case, the difference in energies was compared with the one obtained 
initially. It was found that there is no effect due to varying the 
angle ^NNp and bond length but there is a very significant 



difference when the bond length KNf was minimized. This result shows 
the configuration is affected significantly due to the Rnp stretching 
motion. 





TABLE 13 

Results of experiments for N 2 F 2 (stability) 


Year 

19h7 


Author 
S.H. Bauer 


Reference 

J.Am.Chem.Soc. ,69,310li 


Results 

trans is more stable 


1959 


C.B. Colburn,et,al. 


" •' " " ,81,205 


than Cis configuration 


i960 


R.H, Sanborn 


j.Chem.Phys., 33,1855 


Cis and N=N^p is 
more chemically active 
isomer 


1961 


R, Ettinger,et,al, 


" " •' , 3U,2187 


1963 


G.T, Armstrong 
& S, Marantz 


J.Chem.Phys., 38,169 


Cis is more stable by 
3.0±0.3 ^‘^^Vmole 


Present 






Cis is more stable 
by 3 . 2 U 9 ^‘^^Vmole 



The next interesting calculation is the barrier due to the rotational 
motion which was reported by the Bin^boyn (37) . It was reported by 
experiment to be 32.0 ^°^/mole and calculated to be 8k ^‘^^Vmole (36) 
between trans and linear configuration. However, the present calcu- 
lations were carried out in various ways and those results are listed 
in Table 1$, The calculated barrier between linear and trans form in 
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TABLE 11; 



Equilibriiun Configuration Btn, for N2F2 



Calculation 


Parameters 


APar(^ir -GR ) Barrier” 




Rhh 


R Np (A°) 




Rhu 


Rnf 


(Gir-Ca-ft) 


genera 1 


115 

115 


1.25 

1.25 


l.l;!; 

I.i4i 


0 


0 


0 


0.2553(0, 


II 

Gh 


109.98 

10l;.09 


1.25 

1.25 


l,kh 

1.25 


5.89 


0 


0 


0.2553(0) 


III 

Gk 


109.98 

IOU.O9 


1.2197 

1.2218 


1.25 

1.25 


5.89 


-0.0021 


0 


(-)o.o;^62^) 


nr 

L 


109.98 

IOU.O9 


1.2197 

» 

! 1,2218 


I.2UIO 

1.2l;l;l 


5.89 


-0.0021 


-0.0031 


(-)3,2l;92 

(3.5)_ 


Expt, 

1 value 














(-)3.0i0.3 



Table 15 is too high with respect to the reported barrier but the 
calculation between trans and C2 form shows a satisfactory agreement 
with the calculated (38) and experimental value (39). It must be 
noted that when an attempt was made to calculate the rotational barrier, 
no dihedral angle exists j that is, when the dihedral angle is 90 degrees, 
the bond angle^^-NNfT turns out lineaT, and also the dihedral angle does 
not exist and the motion is predicted to be an up and down motion which 
is shown in Calculation IV and V in Table 15 ♦ The calculated barriers 
are illustrated in Fig. 3. 
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TABLE 1 $ 

Comparison of Barrier for N 2 F 2 



No. Ref. 


Parameters 


Etot. (a.u.) 


ZNNF 


Rnn 


Rnp 


I 


C 2 v (Cis) 


109.98 


1.2197 


I. 2 UO 


-78.87ii0723ii 


II 


C 2 h (trans) 


IOU.O 9 


1.2218 


1.2liiil 


-78.86889607 


III 


0.5(1 + II) 


107.035 


1.22075 


1.2ii255 


-78.8711i8U21 


IV 


02 (N .) 


llU.07 


1.2208 


1.2ii26 


-78.7582U577 


V 


C2;i» ( ) 


180.0 


1.2119 


1.14i0„ 
2£ 


-78.5U333813 



E(II) - E(IV) E(II) - E(IV) E(I) - E(IV) 

barrier ^ 69.U56 = 20U.3^6 = 72.7057 






E(II) — E(V) = 32 kcal/mole (37) 
E(II) — E(V) = 8U kcal/mole ( 36 ) 



Fig. 3 Barrier for N 2 F 2 as function of XNF 
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V-7. Oxygen difluoride (F2O) 

For the F2O molecule with C2V symmetry (RqF — (13)}> the 

lowest energy configuration is calculated to be bent with an/^Op" =97.16° 

(at fixed R^=1,^0A°) as compared with the experimental angle of 103.8° 

(38) . The bond length R^ (at angle of 97.16°) is calculated to be 

1.1792A°. Also the bending and stretching force constants are calcu- 
lated to be l,li;97 X 10“^^ °^S/rad? (observed =0,55 x 10~^^ ®^^/rad^) 

and h0,9h x 10^ '^y^Vcm (observed =5.57 x 10^ ‘^y"®/cm (lU)). 

It is of interest to compare the orbital energies as a function 
of internuclear angle with that publisted by Buenker and Peyerimhoff 

(39) . The calculated results are illustrated in Fig. U and these were 

compared with c alculations by Buenker and Peyerimhoff (using full 
LCAOSCF theory) also in Fig. U. The comparison in Fig. shows a 
satisfactory agreement except for Ib^^ and 3b2 . It is noted that, 

in full LCAOSCF calculation, the lowest doubly degenerate orbitals 
(due to Is electrons of two Fluorines) and the next singly degenerate 
orbital (due to the Is electron of Oxygen) were neglected. Therefore 
the comparison is valid starting with the third level in the Fig given 
by Buenker, 
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V-8. Itydrogen Peroxide (H 2 O 2 ) 



For the hydrogen peroxide molecule, the calcvilations were performed 
to find the barrier and the equilibrium configurations with R^=0.97A°, 
ZH00=100°, R~=1.1 i 9A° and dihedral angle of 93.6° (13)* Unfortunately, 
the present calculation does not show the existence of a dihedral 
angle. The result of the calculation shows that there exists a Cis 
and trans configuration for '^ 2^2 ^ N 2 F 2 . The structiiral 

parameters are illustrated in Fig. 
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V-9. Homonuclear diatomic molecules 



For the diatomic molecules, the equilibrium configurations of H2, 

Liz, Nz and Fz were calculated. An attempt was made to find a new set 

of optimized parameters by changing the Slater exponents and beta 

zeroes (/J ), The equilibrium bond lengths and force constants are 

listed in Table I6, It is noted that Bz, Cz and O2 molecules are excluded 

in the calculations because they are ground state triplets, and the Bez 

molecule ip non-bonded. The necessary partial derivatives were foxmd and 
i 

listed in Table I6. 



TABLE 16 

Equilibrium Configurations and Some Partial derivatives 
for homo nuclear diatomic molecules 





Re (k 


°) 




average 
and aisc 




between (S 


,E)o± 0.1 




Cal, 


Obs. 


Cal. 


Obs . 










H2 


0,7U6 

(0,89) 


0.7U 


10. lU 
(6.5) 




-0.5U 


-0.3it 


19.18 


1.0555 


Liz 


2,179 

(2.789 


2.67 


0.8ii 

(O.hl) 


0,25 


-2.9U 


-0,022 


3.07 


0.52 


Nz 


l.lU 


1.09 


52.16 


23.10 


-0.U7 


-0.01 


61.06 


2.3h 


Fz 


1.12 


iM 


56.U1 


U.U5 











For those molecules, an attempt was made to find a new set of 
optimized parameters fitting experimental equilibrium bond lengths 
and force constants. The Slater exponent and were considered in 

adjusting the parameters to optimize agreement for calculated potential 
surfaces. The procedures used to find a new set of Slater exponents 
and beta zero ( ) were as follows: 

First: Change the Slater exponent by i 0,1 and beta zero by il,0, 

then find the calculated equilibrium bond length and force constant. 
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Second: 



Third: 



Fourth: 



Find the partial derivatives 






^S JyS 









•^R 






dtnA. 









vjhere ”S", stands for Slater exponent. 

Find the new Slater exponent and by following equations: 



dh-. (-||-Vd5 + (^)s’46 

then approximate the dR, ds and 



as follows: 



(5.2) 



C|R = Re^.— 

ds — = Snew — Scisl- 



(5.3) 



=yf3 — j 3 cfi\. 

Solve equation (5*1) and (5.2) to find the new Slater exponent 
and . For example: 

(1) H2 molecule: 



-0.06' 




'-5.1: 


-3.U 




-0.277' 






.-U.83 

/ 




19.18 


1.055 




O.U57 




0 

— 



theref ore 



(2) Li 



Snew = 0,923 
jSnev = ^. 1:57 



(5.U) 



2 case: 



'-0.1:9' 




2.91: 


0.022 




f 

-0.161: 






-0.59 




3.07 

■V 


0.52 

J 




-o.Uio 




0 

1 



therefore 



(3) N2 case: 

-0.05 ' 
-29.06 

therefore 



Snew = 0.1:9 
^new = 8.59 



(5.5) 



-0.1:7 


-0.01 




0.87 




'^5 ' 


61.06 


2. 31: 




-35.5 




0 

\ 



Snew =2.87 
^new = -9.5 



(5.6) 



Fran these results, (5,U)and (5.5) are consistent but (5.6) 
is completely inconsistent because the force constant is too 
high. 

Fifth; Finally, repeat the calculation using the above new set of 
input parameters (Snew and jQ new). The result for the H2 
and Li2 molecules are listed in parentheses in Table I6, 

VI. CONCLUSION 

The Approximate Self-Consistent Field M.O. theory was applied 
to the calculation of electronic structure for several selected mole- 
cules and these results were compared with results by full LCAOSCFT 
calculations, if they were available. The most important features of 
these calculations are to predict equilibriiun configurations (bond 
length and bond angle) and the barrier by the variation of internuclear 
distances and angles using the ASCFT, The EHT calculations are limited 
tc hydrocarbons which give good results (refer to calculation of Ethane 
molecule). In general, the barrier predicted by ASCFT calculation is 
in satisfactory agreement with the experimental value. 

An attempt was made to find a new set of input parameters (Slater 
exponent and ) to improve the prediction of potential surfaces 

for certain homonuclear diatomic molecules but the attempt failed 
except for the H2 and Li^ molecules (see Table I6). 

The features of Potential Surfaces of molec\xles can be predicted 
with the approximate self-consistent field theory. Qualitatively, the 
agreement between calculations and experimental observations is quite 
good. It is conceivable that, with the adjustment of certain parameters 



such as resonance integrals ( ), Slater exponents and penetration 

integrals, quantitative agreement with experimental results can be 



obtained 
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I Computer program for Approximate Self-Consistent Field 
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WRITE (51.101) 

DO 20 I = 1.3 
DO 20 J = l.N 

20 WRITE (51.120) (COo(I.J.i<). k = l.N) 
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^SUBROUTINE INTER ( N »NU,C.P .PA, F ) 

DIMENSION 0(40.40), P(40.40). PA(18), F( 40,40) 
WRITE (SI, 100) 
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DO 30 1 = i,N 
Do 3o J = i »N 

KriO(I»J) = u, p* ( Z5 ( I ) +Z0( J ) ) *raD ( I » J ) 
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CALCULATt PEKIOD Z ATOi^ OVERLAP INTEGRALS 



ISUb = 4*1 
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CA(3*i) = CA(1,1)*CA( 1,3)/DEN 
19 DO 15 K = 1*3 
DO 15 L = 1,3 
15 CC(K,L) = CA(K,L) 
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CALCULATE HYDROGEN OVERLAP INTEGRALS 



DO 60 J = I 1 ,N 
JSUB = NB+J 

60 S(ISUB»JSUB) = SF( 1»RH0( I,J) »T*U{ I ,J) ) 
DO 80 I = 1,NT 
II = I+l 
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60 SF = (P**5/16.0)^;-SQRTF( (1.0-T^:-T)**5) 

1 (ABF(2»2,P,TS)*(ABF( 1»0»P,TS) + ABF ( 1 » 4 , P , TS ) ) 
2ABF( 1,2»P»TS)*(ABF(2,0,P,TS) + ABF ( 2 » 4 , P , TS ) ) ) 

return 

70 SF = (P**5/32.0) * SQRTF( ( 1.0-T*T) **5 ) * 
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I 



12 FACT( I ) = FACT! I-l )*FI 
IFBF = 12345 
20 ABF = 0.0 

GO TO (30»130) KAB 
30 DO 40 M = 1»KP1 
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G(J,I ) = G( I ,J) 

RD = RAD( I » J) 

P = ZS( I )*RD 

V( I,J) = VF(2»P»RD)*2(J) 
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IF(IPT(4)-2) 90,85,90 
85 DO 87 I = 1,N 
DO 87 J = 1,N 
87 V( I ,J) = Z( J)*G( I ,J) 
DO 88 I = 1,N 



FUNCTION GF ( KEY»P»T» R ) 
IF(ABSF(T )-0. 000001 ) 1 ,1,2 
GO TO (110,120,130) KEY 
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10 VF = (1.0-(1,0+P) *EXPF (-2.0»P ) ) /R 
GO TO 100 

20 VF = ( 1.0-( 1.0+1, 5*P+P**2+P**3/3. 0 )*EXPF(-2.0*P) ) /R 
100 RETURN 
END 
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CALCULATES EIGENVALUES AND EIGENVECTORS OF REAL SYMMETRIC MATRIX 0009 
BY THE MODIFIED GIVENS METHOD. THE PARAMETERS ARE... 000 

NX ORDER OF MATRIX OOlO 
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DIMENSION B(NX>5)» A(l)» ROOT ( 1 ) , VECT ( NX »NROOTX ) 0047 

DATA (KILLBIT = 7777777777777776B ) ♦ ( ET A = l.OE-10) 0048 
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805 VECT(I) = VECT(I-l) 4 0.013125901 0262 
C DON-T GET EXCITED. THE CONSTANT IS JUST A FACTOR TO TRY TO 0263 
C DESTROY CYCLES IN THE STARTING VECTORS. THIS HELPS OUT THE 0264 
C ALGORITHM IN SOME CASES WITH MANY DEGENERACIES. 0265 
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890 ELIMl = VECT(1,I) 0^32 

DO 830 J=1»NM1 0333 

TEMP = B(J,5).AND.1B 0334 

IF (TEMP) 840,850 0335 
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DO 910 J=1»K 0368 

910 TEMP = TEMP + VECT(J,n**2 0369 

TEMP = 1 .0/S0RTF( TEMP ) 0370 

371 

DO 920 J=1»N 0372 
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970 TEMP = TEMP + B ( J » 2 ) * VECT ( J » K ) 0403 

TEMP = TEMP + TEMP 0404 

DO 980 J=IM»N 0405 

980 VECT(J,K) = VECT(J»K) - TEMP*B(J»2) 0406 



r- 00 c^ o rH c\j cn 

<t <t <t <t 



o 


c 


O 


c 


• o 


o 














• 


















































K 


















































l-H 






































o 












IL Z 






































(X 












o 






































o 














CO 




































o 












LU Q 


z 






z 






























u 




9 








CO < 


o 






< 










































Z> LU 


< 1- 






l-H 






























z> 




00 








DC 


< 






CO 


































f-H 


9 






cO 


z 






LU 






























< 












LU LU 


•-H lj_ 






1— 






























1— 




_J 


Ch 






X CD 


o 






CL 






























LU 




o 


•w 






< 


CD 






< 






























CD 




z 


CO 






21 X 


LU CL 






U 






























9s 




O 


M 






< 


> LU 




































lO 




CD 


cO 






O X 


O CD 






o 






























rsj 












o 


z z 






z 






























0s 




9 


•v 






Z CO 


LU X 




































rsJ 






■ — . 






U 21 


CL Z 






H- 






























9- 




00 


C^ 






O 








u 






























_J 




t— H 


W- 






1— 1— 


CL LU 






3 






























o 




•W 


NJ 


9 




< < 


LU X 






DC 






























z: 




X 


CO 






X 


1— 1— 






1— 






























o 




CL 




00 




1— >- 


< 






CO 






























CD 






9 


r— 1 




21 


_J o 






Z 


































9 


*— 






X X 


H- 






o 






























X 






00 


h- 




U 3 


LU 






u 


• 




























Q_ 




00 


1 — 1 


< 




^ o 


CD _J 








CO 




























9 






9 


o 




CO 


< 






z 


LU 




























X 






CO 


z: 




X 


-J ZD 






l-H 


U 












cO 
















1— 




X 








1- _J 


_J O 








z 












X 
















«K 




1— 


o 


9 


00 


-• _J 


l-H LU 






Q 


< 












X 
















_J 






DC 




f-H 


< 


IS 






LU 


h- 












< 
















< 




9 


O 


Ch 


9 


cO 


1— 






CO 


cO 












CO Z 






















O 


'W' 


CO 


LU 1- 


>- LU 


• 




X 


l-H 












X -H 
















a 




oo 


u 






o -- 


LU iO 


CD 


o 




o 












< CD 
















1— 




1 — 1 




< 


LU 


z: z: 


X 


LU 


CD 


CO 














l-H CL 
















CO 




W 


9 


CO 


u 


<L 


1- X 


l-H 


z 


1- 


o 








CO 




X O 
















*—{ 




_J 








(X 


-J 


a 


u 


z 


z 








h- 




Z O 
















9 




< 


O 




9 


cx • 


-J 


ZD 




l-H 


< 








X z 




X u 


cO 














h- 








— 


— 


< z: 


• < 


u 


LU 


o 










I— LU 




t— 


t— 














Q_ 










00 


l-H 


o -- 


u 


LU 


CL 


CO 


< 


< 




O Z 


CO O z 


X 














►-H 




— 






f-H 


CD 


1- O 


CO 




LU 


X 


1— 




Z CO O 


O 


X < 


o 


















00 




CM 




Z CD 


II 






o 


_J 


Q_ 


LU 


►—1 


LU LU Q_ 


X 


l-H 


l-H 














o 




f-H 


9 


w 


1— 


< < 


z 


CO 


I 


z 


KD 




X 


X 


_J O X 


X 


Z CO 


X 














z: 








0. 


CO 


LU 


CL *-* 


_J 




l-H 


Z 


< 


1— 


0. 


CL LU 


M 


O X 


:s 














9 




0. 


o 




l-H 


< CL 


LU 


< 


CO 


1— 


< 








ID < 




•-H 1— 
















X 




1— 




CO 




1— 


CD CO 


H- 


4^ 


X 




CO 


CO 


cr> 


Z X X 


< 


X X 


u 














z 




cO 


w 




9 


< CO 


Z =) 


l-H 


o 


< 


o 


LU 


LU 


LU 


O U LU 


h- 


< < 


l-H 














r> 




►— H 


< 


9 


-<— • 


CD 21 


^ X 


CD 


l-H 


1— 


z 


-J 




_J 


CD 1- 


X 


fsj U 


z 














z: 










CO 


O 


Z 1- 


IX 


1— 


CO 


o 


o o 


O 


LU < 


CD 


l-H 


o 


















9 


LU 


Ch 


9 


Z 1- 




O 


CL 




CD 


z 


Z 


z 


Lu X _J 




Z X 


X 




►— < 










1— 




— - 


CD 


9 


CO 


l-H < 


U CO 




C Ll_ 




< 


< 


< 


O O CO 


X o o 


< 














D 




m 




CM 






*-H l-H 


Li- 




o 


z 








U 


o 


l-H 






I 










0. 






9 


— ' 


00 


CO li- 


z 


O 


Lu 




o 


lu 


lu 


tu 


X X 




X 


X 














2 




1— 




CO 


r-H 


CD O 


o z 




O 




CL 


O 


O 


o 


< X O 


X 


X < 


o 




0. 


r— 1 












D. 


00 ZD 




< 


1— 


IX 




< 


Lu 








X O 


< 


O X 






5: 














HH 


r— 1 t/) 


X 


LU QC 


< 


LU 


X 


CL 




X 


X 


X 


X X 


X 


X 


X 




z> 


1 








LU 






— 




O DC LU 


• 


CD 


< 


CL 


CO 


< 


< 


< 


< X < 


X 


X < 


< 


LU 






LU 






Z 




z 


CO 




DC 


CD 


ct: LU z 


CL 


< 


Ul 


CL 


CL 


CL 


< X 


< 


< 


X 


=3 






z> 






►—H 




o 


M 


— ^ 




LU 21 


z 


D 


(L 




H- 


CL 


CL 


CL 


XXX 




X X 


X 


z 


II 




z: 


z: 




I— 




l-H 








z: z> 


LU 


Z 


< 


LU 


< 


< 


< 


< 


X X < 


X 


X X < 


t—t 








DC 




D 




CO 


9 


' — 




z 


X 1- 






X 


Z 








h- < 


X 


< X 




h- 


CL 


II 


1— 


r> 




O 




z: 




< 


00 


1— 


1- ID 


LU 


Z 


h~ 


l-H 


LU 


LU 


LU 


X 


1— 




X 


z: 


51 




z: 


I— 


o 


DC 




LU 


00 


1— 


f-H 


3 LU 


O 


X 


< 




ID 


X 


X 


X 


CO X X 




X CO 


X 


O Z) 


5: 


o 


LU 


z 


X) 


< 


z: 


f-H 


LU 




O X 


X IX 1— 




CO 


CL 


1— 


1— 


1— 


-H X H- 


CO 


X -- 


X 


u 


“) 


» — * 


u 


DC 


LU 






►— 1 


w 


CD IS 


CL \- 


1— CD 




cO 


►— H 


O 








X 


»— 1 


1— 
















J) 




O M 


cO 


< 




ZD 


cO 


•— 1 




O 


i/) 


CO 




-J CO 




O /) 
















f—H 




1— 1 


CM 


CO 


CO CO 


IS CO 


•-H 




Q_ 


U 


HH 


l-H 


l-H 


ID CO -H 


< CO X 


l-H 


























l-H l-H 






1— 


1— 










z 


h~ 


o 










CNJ 


f— H 
















X 




o 


CL' 


cO 




-J 


X 


X o ^ 


X 


o 


IS« 


o 




o 


O 


o 
















z 




z 


•—I 


l-H 




< h- 


CL 


CD NJ Nj 


CO 3 LJ 


<■ 


vO 




m 


o 


o 










































o 




o 


r—H 


1 — 1 


u 














u u 


U U 


u 


U 


u 


U 


c 

c 


u 


c 

c 

c 


U 


c 

c 


u 



85 . 



I 



THE FOLLOWING ARRAYS ARE READ IN ONCE AND ape USED TO CONSTRUCT THE 
ABOVE ARRAYS EACH TIME AN ATOM IS FED IN 
SZ IS ARRAY OF CORE CHARGES 
SZS IS THE ARRAY OF SLATER EXPONENTS 
SBETA IS ARRAY OF BETA ZEROS 
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121 FORMAT (22H TABLE OF CORE CHARGES) 

122 FORMAT (26H TABLE OF SLATER EXPONENTS) 

123 FORMAT (36H TABLE OF BETA ZEROS» NEGATIVE UNITS) 

12^A FORMAT (AlH TABLE OF lONIZ. POT.> S ORBITALS, CNDO/ 1 ) 
125 FORMAT (AlH TABLE OF lONIZ. POT., P ORBITALS, CNDO/ 1 ) 



12b FOKriAT (39rl TAdLc UF U • £)* ( I +A ) , S ORbITALo» CTnIDO/Z) 

12/ FOKhAT (39H TAoLL OF 0.5--(I+A), P ORbITALS» ChDC/2) 

126 FORMAT (23H TAoLE OF ATOMIC MASSES) 

129 FORMAT (67H STP AT. NO. ALPHA THETA PHI 

IBOND LENGTH ) 
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DO ^73 J = 1,IA1 
273 CD( I »J) = CE( I »J) 

DO 273 I = 1»IA1 

275 CD(1»I) = CD(1,I)-B0NDL(1A) 
280 DO 283 I = 1»3 
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44 ) = CD(K»J) 

ISN = NOAT(J) 

JS = J+1 
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ETHANE (ECLIPSED) 
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APPENDIX 



Computer output of Water (H2O) 
(Sample calculation) 
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UISIANCES IN atomic units 

0 1. 9^4377 1.944377 
1.944377 0 3.068970 
1.944377 3.068970 _0 



0 3.37r355 3.370355 
.378355 0 Ai682764 
.378355 3.662764 0 
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ohbital energies 

-1.442662 -.75990^> -.701810 -.654866 .289195 .309628 
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